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Abstract 

We investigate the behaviour of a system of particles with the different 
character of interaction. The approach makes it possible to describe systems 
of interacting particles by statistical methods taking into account a spatial 
nonhomogeneous distribution of particles, i.e. cluster formation. For these 
clusters are evaluated: their size, the number of particles in a cluster, and the 
temperature of phase transition to the cluster state. Three systems are under 
consideration: electrons on the liquid helium surface, particles interacting by 
the shielding Coulomb potential, which are found under the influence of an 
elastic field (e.g. nucleons in a nucleus), and gravitating masses with the 
Hubble expansion. 

Key words: statistical mechanics, interparticle interactions, clusters 

PACS: 34.10+x General theories and models of atomic and molecular col- 
lisions and interactions (including statistical theories, transition state, stochas- 
tic and trajectory models, etc.) - 36.90.f Other special atoms, molecules, 
ions, and clusters. 



1 Introduction 

Basically, statistical description of many-particle systems concerns homogeneous 
states. In monograph [1] , a general method based on quantum field theory was suc- 
cessfully applied to homogeneous states of condensed media. An approach that took 
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into account spatial nonhomogeneous states of a system of particles was first pro- 
posed in paper [2]; the approach was resting on the employment of the methodology 
of quantum field theory as well. The most general approach to the description of a 
system of particles, which included the interparticle interaction was offered in paper 
[3] in which the procedure of discount of the number of states providing for the 
nonhomogeneous particle distribution was proposed. However, a correct solution of 
nonlinear equations in the saddle-point is an impossible task for the most part: the 
general approach does not permit one obtains analytical solutions for a system of 
particles when the inverse operator of the interaction energy does not exist. In the 
interim we need to know a series of parameters which characterize possible clusters. 
The major parameters are the cluster size and the critical temperature of cluster 
stability. Because of that, one can arise the question about a more simplified ap- 
proach, which, however, would allow the study of cluster formation in an initially 
homogeneous system of particles. 

In the present work we study a system of interacting particles having regard 
to the concept proposed in papers [2,4]. Such an approach reduces the number 
of variables, which determine the free energy and introduces a peculiar "combined 
variable" [5] playing a role of a typical variable that describes the thermodynamics 
of the system with a spatial nonhomogeneous particle order. 

For instance, as an analog of a certain combined variable one can consider the 
interaction of two waves on the water surface, which interfere in a limited range. 
However, after the interference the waves depart and hence in this case the combined 
variable does not bring into existence any instability of the system studied. On 
the other hand, the interaction of an incident electromagnetic wave with a polar 
crystal results into the polariton, a stable formation, which is indeed described by a 
combined (canonical) variable - a mixture of the electromagnetic field and optical 
phonon. 

Can the combined variable give rise to the change in the homogeneous distri- 
bution of interacting particles? We analyze such an option considering concrete 
examples. Below we take a look at the specific systems with long-range repul- 
sion/attraction and short-range attraction/repulsion. In particular, the approach 
provides a way of estimating of the cluster size, the number of particles in a cluster 
and the temperature of phase transition to the cluster state. The approach is exam- 
ined in three examples: i) electrons on the liquid helium surface, ii) particles with 
the shielding Coulomb interaction, which are found in an applied elastic potential, 
and iii) gravitating masses with the Hubble expansion. 

2 Statistical mechanics for model systems with 
interaction. The principle of states selection 

Let us consider a model system of particles with arbitrary kinds of interactions. 
Nevertheless, the de Broglie thermal wavelength A of a particle is taken to be greater 
than the mean distant I between particles. However on the other hand, the value of 
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A must be smaller than the mean amplitude of scattering. The type of the statistics 
introduces significant peculiarities to the interparticle interaction, which, however, 
can be rendered by the classical methods allowing one does away with dynamic 
quantum correlations. In this case the macroscopic state of the system in question 
may be specified by filling numbers Us, where the symbol s describes the state of 
the system. In other words, particles can occupy only knots are still non-occupied. 
This means that the filling numbers Ug can run only two meanings: 1 (the sth knot 
is occupied) or (the sth knot is not occupied). Besides, the allocation of particles 
in the lattice invites the type of statistics, which set limits on the behaviour of 
particles, i.e. the statistics refers particles to either bosons or fermions, see below. 
The Hamiltonian of the system of interacting particles can be written in the form 

'^("') = ^^'^^ Vss'nsUs' + ^ Uss'TlsUs' (1) 

s ss' ^ ss' 

where Es is the additive part of the particle energy in the sth state; Vss' and Uss' are, 
respectively, the paired energy of attraction and repulsion between particles in the 
states s and s'. Signs before the potentials in expression (1) reflect the proper signs 
of the attractive and repulsive paired energies and therefore both of the functions 
Vss' and Uss' are positive. 

Such writing of the Hamiltonian corresponds to the model in which particles are 
found in sites of the three-dimensional lattice. As will be evident from the further 
consideration, a specific form of the lattice will be unneeded. However the medium 
will be supposed isotropic at the transition to the continual approximation. 

Note that the Hamiltonian of the model of replacement of solid solutions [6] 
amounts to the Hamiltonian (1) as well. This means that the approach proposing 
can describe the behaviour of defects in solids, allows, liquids, etc. (in these cases one 
particle will characterize an unfilled place and another particle that of a filled one). 
The reduced solution Hamiltonian can be amounted to a general Ising model, which 
takes into account an arbitrary interaction allowing the spatial-nonhomogeneous 
distribution of particles. 

The partition function of the system is 

Z^Yl eM-'H(ns)/kBT) (2) 

{Us} 

where the summation provides for all possible values Us, i.e., all states of the system. 
Such summation can be formally performed, if one passes to field variables according 
to Ref . [2] . It is known from the theory of Gauss integrals that 

jj^ r \ 

*^^P[y Z] 'UJss'risns'] = Rc y Dx exp[z/ ^ '^sXs " 2 ^ '^ss^sXs'] (3) 

ss' s ss' 

where Dx = H^ ^ det 1 1 Wss' \ \ V^dxs implies the functional integration with respect 
to the field x ^-iid — il in relation to the sign of interaction (-1-1 for attractive 
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and —1 for repulsion). Besides the obligatory condition 

W~],Ws>:s' = Sss' (4) 

s" 

is full-filled. In particular, condition (4) in principle could allow one to determine 
the inverse operator w~) (in the case when it is exist) which in its turn would permit 
the construction of the Green function for interacting particles. 

Let us introduce the dimensionless energies Wss' = Wss'/k-QT and Eg = Es/k-oT. 
Note that hereinafter the tilde over a symbol means the division of the symbol by 

So the partition function (2) can be represented in the form 

oo oo 

Z = He J D(f) J Dip J2 '^[^Y^ ^jig + + i(t>s)ns 

-oo -oo {"s} « « 

-\ J2iKs'^s(l)s' + V~}^sll^s')]- (5) 
ss' 

We do not include the dependence on momentum of particles to the functional 
integrals in the partition function (5); the sum takes into consideration only the 
spatial distribution of particles and their energy. 

Now we can settle the quantity of particles in the system, J2s ^s = N . The 
procedure is an equivalent to the consideration of the canonical ensemble. For this 
purpose one can use the formula known in the theory of complex variable 

— <Ldz z^'-^-^s^^ ^1. (6) 
It makes possible to introduce the sum of the canonical ensemble 

ss' 

+ (7V-l)ln^} X ^exp{^n,(V^, + -In^)}. (7) 

{"4 « 

Summing over Us (note {ug} = 0, 1) we get 

Z = Re^ j D(t) j D^pjdz e^^"^'^'^) (8) 



where 



S = E{-^EfcVA' + *;.'^^^^')+^ln|l + V^^e^«cos0,|} 



s' 



+(A^-l)ln^. (9) 

Here, the function rj is equal to ±1 (Fermi or Bose statistics), see e.g. Refs. [7- 
9]. This provides the means determining the necessary procedure of the states 
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selection that introduces the most essential contribution to the partition function 
and defines the free energy of the system. Let us set 2; = ^ + i( and consider the 
action S on the transit path passing through the saddle-point with a fixed imaginable 
variable Im ^ = (q. In this case, it stands to reason that the action S, similarly to 
quantum field theory, is the variational functional that depends on three variables: 
the fields 0s and t/jg, and the fugacity ^ (here ^ = Q-ti/ksT ^^ere /i is the chemical 
potential). The extremum of the functional must be realized at solutions of the 
equations SS/64>s = 0, SS/Sips — 0, and SS/S^ — 0. These equations appear as the 
following 

e--^-e^-sin(/), 

l^u^As'^-— — i^- -; (10) 

g, ^ + ^ cos 0s 

g, C + V^ -^e^^ cos 0s 

E (12) 

g, ^ + rje ^a'e^s' coscpgi 

Equations from (10) to (12) completely solve the problem of the statistical descrip- 
tion of systems with any type of interaction. One or another of states of the system 
is realized in accordance with the solution of the nonlinear equations. Among the 
solutions there are solutions which correspond to the spatial nonhomogeneous dis- 
tribution of particles. So the Bose condensation is realized in the real space because 
we can say that an aggregation of particles is the Bose condensation of a sort. One 
of the possibilities of cluster formation was demonstrated in Ref . [2] ; the nonhomo- 
geneous distribution of particles which were characterized by the shielded Coulomb 
interaction for both the attraction and the repulsion just yielded stable clusters. 
Such type of the interaction was chosen owing to the existence of the inverse oper- 
ator of the interaction energy in an explicit form. Nevertheless, it is a very difficult 
problem to find the inverse operator among the most widespread types of interac- 
tions. This is why, there is the necessity to develop methods which are capable to 
describe the particle distribution at any interaction. In Refs. [4,5] the approach 
based on the selection of states in the form of a "combined variable" has been pro- 
posed. The combined variable has given us the chance to find connections between 
different fields (such as 0s and '0s in the present paper) along the extremal path 
that passes through the saddle-point. 



3 Combined variable and the selection of states 

Let us introduce the designation 

e^'^^e'^^ COS0S , , 

Ls = = — — . (13) 

{ -I- rje^^^e'^^ cos 0s 

The value Fs includes the two field variables, ipg and 0s, and the fugacity ^. Inserting 
Fs into the left hand side of Eq. (12) we will see that the total sum J2s is equal 
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to the number of particles in the system studied 

^r, = iv-i. (14) 

s 

This means that Tg directly specifies the quantity of particles in the sth state. 
Hence the value Tg might be considered as a characteristic of the number of particles 
contained in the sth cluster. So it is quite possible to consider Tg as the variable of 
particle number and, on the other hand, it may be called the combined variable owing 
to the fact that it combines the variables ipg, 4>s, and ^. Probably the introduction 
of the combined variable is a rough approximation, but it admits to determine 
the clusterization conditions in cases when the inverse operator of the interaction 
potential can not be found. 

Utilizing the variable (13) we can represent the action (9) as a function of only 
two variables, Tg and ^, and it makes possible to do away with the inverse operators 
VgJ: and ujg, in the final expression. 

If we multiply two hand sides of Eq. (11) by the function Vgg" and then summing 
Eq. (11) over s, we obtain 

= X^^'^s'rs- (15) 

s 

Now let us multiply the same equation (11) by ■^s; then summing it over s we acquire 

J2Ks^^s'i^s = ^Vgg^Tg^Tg. (16) 

ss' ss' 

Multiplying Eq. (10) by UgHg and summing it over s we have 

2^ Ugg,Ug''g(l)g' = - 2^ Ugg,> — ^— . ( 17) 

gg, s i + "^^'^ COS (pg 

Prom relationship (13) one obtains 

^e-^=e^= COS0, = . (18) 

Let US substitute (18) into (17) taking into account Eq. (15) and using condition 
(4) for the left-hand side of Eq. (17). We gain 

<Pg = -J2us'sle-^^e^^"'^^"^^^"{l-riTg,)sm^g,. (19) 

s' ^ 

With the formula sin^^/ = ~ cos^ 0^/, cos 0s can be substituted from (18) into 
(19). As a result, instead of (19) we have 



0. = -J2ugJ^e-^^^'e'^^"'^"^^^"{l - riTg,y - T^. (20) 



Consequently using Eqs. (10) and (13) we get 

^ u~g,(j)g>(i)s = - ^ Ts(j)s tan (j)g, 
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or in the explicit form 



ss' ss' ' ^ 



]_g-2£,g2^^„ii,,„r,„^^ _ ^r^)2 _ r2 . (21) 



here, allowance is made for the negative sign of the field 0^ (see (19)). 

Now using relationships (16) and (21) we can rewrite the action (9) in the point 
of extremum as follows: 

- lj2^ss'rs'rs-vJ2^T^\'i--V^s\ + {N-l)ln^-lj2^ss' 



X ^[le-2^»e2E."*.."r."(i _ ^r,)2 - Tf]. (22) 

It has been shown in Ref. [5] that the minimum of the free energy based on expression 
(22) is lower than that obtained in the framework of the mean field approximation. 
Indeed, let us assume following Ref. [7] that Nu and are numbers of particles 
getting into the influence of the repulsive and attractive force respectively. The 
number of particles Fg in a cluster can be put constant equals F. Let interaction 
energies in the range of influence of these forces be change to their average values. 
In this case expression (22) becomes 

S ^ --{N - l){N^rur + N^TVT + 2r]ln\l ^ r]r\} 

--{N - 1)[N^U{1 - 77r)e-2^e'^''^^'r' - Sln^] (23) 

where U and V are average values of the repulsive and attractive potential respec- 
tively. Here, in the mean field approximation the expression in the curly brackets 
coincides with the free energy. The next term, i.e. the expression in the square 
brackets decreases the free energy at ^ < 1. This inequality, ^ < 1, corresponds 
to the coherent state of a system of interacting particles. Thus the approach that 
is developing can be considered as a general methodology of the mean field theory 
applying to the nonhomogeneous particle distribution. 

Expression (22) is a function of only one variable F and the fugacity ^ and it 
is applicable to any kind of the interaction, even though the inverse operator is 
unknown. The most advantage of expression (22) lies in the fact that it provides a 
way of finding of such characteristics of the system as its size, the number of particles 
in a cluster and the temperature of phase transition. 

Now let us pass to continual variables into the action (22). Inasmuch as we 
are interested in the nonhomogeneous distribution of particles in an indeterminate 
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volume, let a radius R be the fitting parameter of the system studied. Assume that 
the density of particles is distinguished from zero only in the cluster volume. Then 
Eq. (14) can be written as 

TK = N-1 (24) 

where K is the number of clusters in the system and F defined in expression (13) is 
the combined variable of the fields (ps and ips, and the fugacity ^ in a cluster. As is 
seen from expression (24) the variable F can be interpreted as the mean quantity of 
particles in a cluster. Then the passage to the continual presentation is realized by 
the substitution 

T.fs-K^ I f{f)df; (25) 

cluster 

here the integration extends for the volume of a cluster and V = is the effective 
volume occupied by one particle where g is the distance between particles, i.e., the 
lattice constant. 

Thus in the continual presentation and with allowing for the assumption (14), 
we can transform the action (22) to the form 



^— j df J dr'v{f — r')T{r')T{f) 



^ j df j df'u{f- r')V{f')V{f) 



-2E(r)p2 / t;(r-r")r(r")dr" J_/^^ \2 



+r}^J\n\l-r]T{f)\df+{N-l)lni (26) 



where the integrals are the same as in Eq. (25). 
The integration is effected by the rule 



- J f{f)df = 4^ J j ^^^^d9 j f{r)r'^di 

cluster 3 5' Q g 



1 ^ 

— injfirydr; (27) 
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We shift here the hmits of integration from and R to g and R respectively as it 
allows one to eliminate the singularity of the integrand (the same procedure is made 
below in expression (29)). Thus, we can now normalize the integrals to the number 
of particles F in a cluster, that is. 

Ft 

i / df^ , 47r fr^dr ~J ^F-l^F; (28) 

VJ (47r/3)^3 J g3 ^ \ J 
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the last approximation is in fact correct as we assume that F satisfies the inequahty 

r> 1. 

Relation (28) allow us to introduce the dimensionless variable x = r / g m integral 
(27). Thereby the rule of transformation from summation to integration, i.e. (25), 
becomes 

pl/3 

;^E/^ 3 j f{gx)x'dx (29) 
^ 1 

(once again, in expressions (27) we shift the limits of integration from and Rtog 
and R respectively as it allows one to eliminate the singularity of the integrand; the 
same shift is made in expression (29)). 

Having integrated the action (26), we should exploit the following relationships 



— J df J dr' v{f — r')r{f)r{r') = — J drv{r)T{f) ■ j dr'T{r') 

2 ^^^^ 

— j dr v{r)T ■ T j dr' — — ■ j dr v{r) = SF^ j dx x^v{gx); (30) 



T/2 



1 

pl/3 



— J df J dr'' u{f — r'')r{r)T{r'') = — j df u{f)T^ = 2>T^ j dxx^u{gx); (31) 

1 

With the transformations, we have used the step function: r(r) — r'd(r — r^ / g^) 
where ^{V - r^/g^) = if r^/g^ > T and i9(r - r^ / g^) = 1 if r^/g^ < V. Besides, 

^ j dfj df'u{f, f')f{f)fir') ^ ^1 dfuif)fir) . J dr'fir) 

pl/S 

= Sf j dxx^u{gx), (32) 



where / = ^Jexp[-2 E + 2r J dfv(f}]. 
Now if we introduce the designations 



pi/s ri/3 



a — 3 J dxx^u{gx); b — 3 J dxx^v{gx), (33) 
1 1 

we will represent the action (26) in the following final form 

S^K.{^-ia-b)r'-~{l- r?r) V^™ + 77ln |1 - rjr\} + {N - l)\n^. (34) 

If we minimize the action (34) by the variable F, we will be able to define the number 
of particles in a cluster, the cluster size, and the temperature of phase transition that 
triggers the nonhomogeneous distribution of particles. It must be emphasized that 
the definition of the parameters of nonhomogeneous formations requires only explicit 
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forms of potentials of interparticle interactions. Then the system itself will select the 
realization which will provide for the minimum of the free energy. Nonetheless, we 
assume by first the availability of clusters and then the conditions and parameters 
of their existence are defined. 

Studying the phase transition of clusters to the spatial homogeneous distribution 
of particles, we should take into account the behaviour of the fugacity ^ = g-M/fesT' 
Indeed, in the classical system the chemical potential 

fx = kBT\n{\%n) (35) 

where n is the particle concentration and the de Broglie thermal wavelength of a 
particle with a mass m 

At = h/^3mkBT. (36) 
Thus in the classical case the fugacity 

^ = A^n « 1. (37) 

In the case of a quantum system the strong inequality docs not hold; in this case 
^ < 1 and ^ is only a slight less than the unit. Thereby when we investigate the 
action (34) searching for clusters at the temperature T < Tc, we may neglect the last 
term as | ln^| is rather small in comparison with terms which inchide highest orders 
of r. However, if we consider the action (34) scanning for the critical temperature 
T = Tc, we should retain the last term because in this case ^ trends to a very 
small magnitude such that |ln^| >> 1 and yet other terms decrease owing to the 
diminution of the value of F. 



4 Clusters in condensed media 

In this section we apply the methodology developed above to systems of particles 
with long-range repulsion (attraction) and short range attraction (repulsion). 

4.1 Electrons on liquid helium surface 

The energy of electrostatic repulsion between electrons can be written as Ugs' = 
l/Aireo ■ Q^/{rs — t^/) where Q is the effective charge of the electron in the helium 
film. The attraction between electrons caused by the deformation of the helium 
film can be taken as Vgs' — |7fiim(?^s — i^s')^ where the coefficient 7fiiin characterizes 
the elasticity of the helium film in respect to the deformation caused by imbedded 
electrons (see, e.g. review [10]). 

In this case the parameters a and b (33) of the action (34) are transformed to 
the following 

ksT J X Sneo gk^T ' ^ ' 
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pl/3 



or 



where 



a = ar2/3; 6 = ^r5/3 (40) 

and a » f3 (hereafter the tilde over a symbol means the division by k^T). Since 
for electrons rj = +1 and the inequality E = E/k^T ^ 1 is legitimate (i.e., 
exp{—2E/kBT) <^ 1 where E is the kinetic energy of an electron), the action (34) 
can be rewritten as 

S ^ K ■ [^aV^/^ - ^/Sr^^/^^lnr + rin^}. (42) 

The extremum of the free energy is achieved at the solution of the equation >S"(r) ~ 
0, or in the explicit form: 

^arV3_^;3r«/^+l+l„^ = 0; (43) 

Retaining leading terms in Eq. (43), i.e. highest powers to F, the equation is reduced 
to 

laT''/^ _ H^r8/3 ^ (44) 
and hence the solution to it is equal to 

^^£a^40^ (45) 
11(3 33 eo7fiim 

where n = {^g^)~^ is the concentration of electrons in a cluster. The sufficiency 
is also satisfied, i.e., ^"'(r) > with the solution (45) and therefore it actually 
determines the minimum of the action (42). 

The critical temperature is defined from the condition F = 2, i.e. when the 
temperature is so high that no more than two particles are able to interconnect. In 
this case Eq. (43) turns into 

_4q,.2«/3 + ^/3 -211/3 

knTc — ; — r 

31n{ 

In the first approximation we may allow that in the critical point the fugacity ^ = 
g-/i/A:BT Yeaches its classical meaning (see expressions (35) and (36)). Expanding 
function In ^ in a Taylor series, we obtain In A|,n ~ 3(Arw^/^ — 1)- Substituting the 
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expansion in Eq. (46) and taking into account expression (36) for At, we arrive at 
the following equation for Tc 

, ^ 4 /27rN 1/3 e^-n}/^ 

kBTc^-(—) . (47) 

Now, let assign the concrete numerical values to parameters. Let the charge Q 
be the elementary charge e. Put n — 2.4 x 10^'^ m^'^ and 7fiim = 8.4 x 10""^^ N/m. 
Then setting these values into expression (45) we get the number of electrons in a 
cluster: T fti 10®. For the radius of the cluster we have R = (3r/47rn)^/^ a; 10~^ 
m. The critical temperature Tc estimated from expression (47) is of the order of 100 
K that is far beyond the critical temperature (4 K) of the liquid helium film. Thus 
the thermodynamic condition of the liquid helium stability sets limits on the real 
temperature T^. of the cluster existence. 

The behaviour of electrons in the liquid helium film has been studied by many 
researchers [10]. An experimental appraisal of the mean number of electrons in a 
near-surface bubblon was approximately equal to 10®. The radius of the bubblon 
was estimated as 10~^ cm. In such a manner our qualitative evaluation of the values 
r and R is in agreement with the experimental data. 

4.2 Shielding Coulomb potential 

Let particles be repulsed by the shielding Coulomb potential Ugs' = l/47reo x 
Q2g-K|rs-r^;| y. (^^^ _ where K, is the effective radius of screening of nucleons in 

a nucleus, and attracted by the potential Vss' = ^lii^s — fs'Y that can be created by 
an outside reason (so 7 is the force constant of an off-site elastic field) . In this case 
the parameters a and b (33) of the action (34) become 



pl/3 



4:Ti€ogkBT 



J X 
1 

[ngV^'^ + l)e-^^^"' + («:^ + 1)6-"^]. (48) 



We will consider weight nuclei (F >> 1) and therefore V^^^ should be several times 
as large as the unit. Therefore, at the first approximation we may neglect terms 
with e~'^^^^^^ and then expression (48) becomes 

^^^^2^r^^^9 + l)e-^'. (49) 



Now 



ft ^ 3 4^ / ^2 2^ = - 1^ F^/^ (50) 

2kBT J 10 k^T ^ ' 

Let us rewrite these expressions in the form 

a = a; 6 = ^V^l^ (51) 
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where 

^-J^^^(^9 + l)e-'; P-^^ (52) 
47reo n^g-^k-eT 10 k-^T 

and a » j3. Let us now substitute expressions (49) and (50) for a and h in the 
action (34). With regard for inequahty E » 1 and with allowance for r/ = +1 we 
will get 

S^K-[]^aT^ - ^/3r^^/3 + lnr + rin^}. (53) 

The extremum of the action S is reached with resolving the equation S'{T) = 0, 
or explicitly 

ar-y/3r«/3 + ^ + lne = 0. (54) 
Retaining two highest order terms in Eq. (54), we obtain 



6 a\3/5 



Ul 13^ 



3 *- 11 eoAC 7 



3 \l/3 K ,11 3/5 



^47r^ n 



3/3x1/3 K 



5 ^47r^ 



The critical temperature Tc of the cluster destruction is defined from Eq. (54) if 
one puts r = 2: 

^BTe«^^. (56) 

However, having found the fugacity ^, we now cannot follow the procedure de- 
scribed in the previous subsection. In a nucleus, the temperature and the thermal 
de Broglie wavelength are others. In other words, we should exploit the laws of high 
energy physics, i.e. that the total energy of a nucleon is linked with the momentum 
of the nucleon by the relation E = i/p^c^ + m^c^. Hence in this case 

-kBT = c^p^ + m^c^ (57) 

and therefore 

V ^ 'ikTiT/2c. (58) 

Thus the de Broglie thermal wavelength of the nucleon in a nucleus is determined 
as 

ScksT 

Expanding function In,^ = ln(A|.n) in Eq. (56) in a Taylor series we obtain In^ ~ —3 
and, therefore, equation (53) becomes 

4 

AnQ"^ r / 3 \ i/3 K 1 r / 3 \ i/3 ft; 1 
13 



The results obtained in this section may account for the reasons of the atomic 
nucleus stability from the microscopic viewpoint. The shielding Coulomb potential 
u = l/Aneo x Q'^e^'^''' jr is the typical nuclear (Yukawa) potential that provides for 
repulsion between protons. The potential v = ^'jr^, which is applied to nucleons, 
is ensured their mutual attraction. For example, setting 7 = 4 x lO-*^^ N/m we 
obtain from expression (55) F ^ 30, which corresponds to the number of nucleons 
in a nucleus of zinc. The critical temperature of fission of a proton-neutron pair 
calculated from expression (60) results in ~ 10^^ K. Note that this value is 
several times of magnitude greater than the typical temperature of a weight nucleus 
(the nucleus temperature is of the order of the Coulomb repulsion between protons 
in the nucleus). 



4.3 Gravitating masses with Hubble expansion 

Gravitational physics may also be assigned to the condensed matter [11]. The 
gravitational attraction between particles, i.e. big masses, is beyond question. How- 
ever the uniform character of the Hubble expansion is in doubt. We will use only 
the fact that an additional kinetic energy of particles Eggi = ^m{ws — Wg')'^ is asso- 
ciated with such an expansion. The relative velocity Wg — Wg' of particles which are 
found in points s and s', respectively, pertains to the relative distance between the 
particles, since 

Ws - Ws' = H{rs - Ts') 

where H is the Hubble constant. In this connection, we will consider a model system 
of identical gravitating masses, i.e. stars, with the attraction Vsgi = Grn^ /{vg — Tgi) 
and the effective repulsion Ugg' = \H'^m{rs — Ts'Y- For such kinds of interactions 
one has for the parameters a and h (33): 

pl/3 

""^^^k^ J ""'^"-lo^^^ ' ^''^ 

b-S^7-x^dx^'-%ry^ (62) 

gkBT J X 2 gksT ^ ' 

or 

a = aF^/^; h = ^V"^'^ (63) 

. 3 g^H^m ~ 3 Gm' 

Oi — , ^ ; p = — — . 64 

10 kBT ' ^ 2gkBT ^ ' 

Inasmuch as gravitating masses are associated with the Bose statistics, i.e. rj — —1, 
the action (34) is rewritten as follows 

S K ■ {-ar^^/^ - -/3F^/=^-lnF + Flne}. (65) 



where 
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In expression (65) we omit the exponent term because of the assumption that the 
temperature T of the universe is very close to the absolute zero (at small T the 
exponent in expression (34) tends to zero). Then the equation S'{T) — becomes 

^ar«/^-^^rV3_ i+ine = o. (eg) 

Here, a « (5. Retaining highest order terms in equation (66) we immediately 
obtain the solution 

~ 11 a 33 //2 ■ y"") 

where n— {^g^)~^ is the concentration of particles, i.e. stars. 

The equation for the critical temperature Tc that determines the break-up of two 
coupled masses is obtained from Eq. (66) in which we set F = 2 (and the exponent 
term is neglected). With the inequality a « P, we get 

, _ 22/3 f3 /27rxi/3 Gm^nVs 

By definition, the fugacity ^ = . At the transition from the Bose statistics 

to the classical one, the chemical potential gradually decreases with temperature 
starting from a value // ~ (but still // < 0) to // << 0. So, the function In^ 
in Eq. (68) is strongly negative: In^ = —\yL\/k-QT. This signifies that the critical 
temperature in the system of gravitating masses is absent. In other words, in the 
universe, masses trends to the aggregation at any conditions. 

Thus the results obtained demonstrate that the availabihty of the Hubble ex- 
pansion is the deciding factor in the formation of galaxies in the universe (note that 
without H masses in the universe will trend merely to the total mergence and any 
spatial distribution will not be observed). Actually, expression (67) shows that Rub- 
ble's constant, H 50 km/ (s-Mpsc) = 1.6 x 10~^^ s~\ is one of the main parameters 
that determines the steady state of masses in a cluster, i.e., galaxy. So the greater 
the mean mass m, the greater the number F of stars in a galaxy. 



5 Conclusion 

In the present paper the statistical description of a system of interacting parti- 
cles, which allows the study of their spatial nonhomogeneous distribution, has been 
proposed. One of the most important properties of our model is the presentation of 
the total potential of interacting particles in the form of the paired energy of attrac- 
tion and repulsion. We have started from the Hamiltonian that describes interacting 
particles located in knots of the three-dimensional lattice and applied the quantum 
field theory methodology based on the functional integrals for the description of the 
partition function of systems studied. 

We have shown that the solutions of equations characterizing the behaviour of the 
field variables yield the full resolving the initial statistical task. The solutions may 
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correspond to both the homogeneous and nonhomogeneous distributions of particles 
and only the character of interparticle interaction and the temperature define the 
kind of the solution. 

The free energy represented through field variables in the saddle-point has al- 
lowed us to investigate special conditions, which provide for the formation of clusters 
in a system with the initially homogeneous particle distribution. The conditions are 
answerable also for the size of clusters and the temperature of phase transition to 
the cluster state. 

The approach proposed does not need fitting parameters: the cluster size and the 
temperature of phase transition to the cluster state are completely determined by 
the value and character of paired potentials in the systems of interacting particles. 
Besides, the introduction of a peculiar type of the combined variable has simplified 
the investigation procedure. Specifically, the procedure has provided a way for the 
obtaining the necessity criterion of the transition of a system of separate particles 
into the spatial nonhomogeneous state. 

The authors are very thankful to the referees for the concrete constructive re- 
marks that allowed a significant improvements of the contents of the paper. 
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